DYNAMICAL INVARIANTS AND BERRY'S PHASE FOR GENERALIZED 

DRIVEN HARMONIC OSCILLATORS 



BARBARA SANBORN, SERGEI K. SUSLOV, AND LUC VINET 

' Abstract. We present quadratic dynamical invariant and evaluate Berry's phase for the time- 

, dependent Schrodinger equation with the most general variable quadratic Hamiltonian. 

< 

^ ■ 1. Introduction 



In the previous Letter [23], the exact wave functions for generahzed (driven) harmonic oscillators 
[2], [3], [20], [26], [29], [53], [55], [56] have been constructed in terms of Hermite polynomials by 
transforming the time-dependent Schrodinger equation into an autonomous form [57] . Relationships 
with certain Ermakov and Riccati-type systems have been investigated. A goal of this Letter is to 
find the corresponding dynamical invariants and to evaluate Berry's phase [1], [2], [IQ], [52] for 
quantum systems with general variable quadratic Hamiltonians as an extension of the works [3], 
[12], [IE], [21], [22], [26], [36], [35], [36] (see also references therein). 



Oh. 



'nJ" . 2. Generalized Driven Harmonic Oscillators 

ly^ ■ We consider the one- dimensional time-dependent Schrodinger equation 

o : = (2.1) 

where the variable Hamiltonian H = Q {p, x) is an arbitrary quadratic of two operators p = —id/dx 
and X, namely, 

X ' = -a (t) i',, + b (t) x^ij - ic (t) x^, -%d{t)^-f (t) xi) + %g (t) (2.2) 

(a, b, c, d, f and g are suitable real- valued functions of time only). We shall refer to these quantum 
systems as the generalized (driven) harmonic oscillators. A general approach and known elementary 
solutions can be found in Refs. [3], [5], [6], [7], [lO], [13], [15], [IS], [23], [21], [30], [33], [32], [53] and 
[56]. In addition, a case related to Airy functions is discussed in [25] and Ref. [8] deals with another 
special case of transcendental solutions. 

In this Letter, we shall use the following result established in |24] . 

Lemma 1. The substitution 

i[a{t)x^+5(t)x+K(t)) 

'4' = r=fT ^(^'^)' ^ = /3{t)x + e{t), r = 7(t) (2.3) 
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transforms the non- autonomous and inhomogeneous Schrddinger equation Ii2.2\) into the autonomous 
form 

- iXr = -X^c + cof X (co = 0, 1) (2.4) 

provided that 

— + b + 2ca + Aaa^ = (2.5) 
dB 

^ + (c + 4aa)/3 = 0, (2.6) 



dt 

+ a/32 = (2.7) 



and 



Here 



dt 
dS 

— + (c + Aaa) 6 = f + 2ga + 2coa/3^£, (2.8) 
dF 

- = (^? - 2a6) P, (2.9) 

dK 

^ = gS-a6^ + Coal3h^. (2.10) 
dt 

I u' d . . N 

a = —— . (2.11) 

4a n 2a 



The substitution (12.111) reduces the inhomogeneous equation (12. 5p to the second order ordinary 
differential equation 

/i" - r (t) fi' + Aa {t) /i = Co {2af /^V, (2-12) 
that has the famihar time-varying coefficients 

r(t) = - -2c + 4:d, (T(t) =ab-cd + d'^ + -(--%] . (2.13) 
a 2 \ a d J 

When Co = 0, equation (12. 5p is called the Riccati nonlinear differential equation [50], [51] and the 
system (I2.5p - (l2.10p shall be referred to as a Riccati-type system. (Similar terminology is used in 
[H], [15] for the corresponding parabolic equation.) If cq = 1, equation (12.120 can be reduced to a 
generalized version of the Ermakov nonlinear differential equation (see, for example, [6], [13], [27] . 
[16] and references therein regarding Ermakov's equation) and we shall refer to the corresponding 
system (12.50 - 02. lOp with cq 7^ as an Ermakov-type system. Throughout this Letter, we use the 
notations from Ref. [24] where a more detailed bibliography on the quadratic systems can be found. 



Using standard oscillator wave functions for equation (12.40 when cq = 1 (for example, [T7] . 
and/or [3l]) results in the solution 

i(aa;^+(5x+K)+i(2n+l)7 

ij^ (x, t) = e-(^^+^)'/2 ^ ^) ^ (2.14) 

v/2"n!/iv^ 

where Hn (x) are the Hermite polynomials [39j and the general real-valued solution of the Ermakov- 
type system fl23D- fl230D is available in Ref. [21] — Lemma 3, Eqs. (42)-(48). 

The Green function of generalized harmonic oscillators has been constructed in Ref. [1]. (See also 
important previous works [11], [31], [53], [56], [57] and references therein for more details.) 
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The corresponding Cauchy initial value problem can be solved (formally) by the superposition 
principle: 

/oo 
G'(x,y,t)^(y,0) dy (2.15) 
-oo 

for some suitable initial data (x, 0) = (x) (see Refs. [1], [12] and [16] for further details). The 
corresponding eigenfunction expansion can be written in terms of the wave functions fl2.14p as 
follows 

oo 

(x, t) = ^ c„ (x, t) , (2.16) 
where the time-independent coefficients are given by 

ip*^ (x, t) ip {x, 0) dx 



2 • ^^-^^^ 

\ip^{x,0)\^ dx 

-X) 

This expansion complements the integral form of solution (12.151) . 

The maximum symmetry group of the autonomous Schrodinger equation f l2.4p is studied in [37] 
and (see also |19] and references therein). 



3. Dynamical Invariants for Generalized Driven Harmonic Oscillators 

A concept of dynamical invariants for generalized harmonic oscillators has been recently revisited 
in Refs. [6] and [16] (see [9], [10], [11], [31], [32] and references therein for classical works). In this 
Letter, we would like to point out a simple extension of the quadratic dynamical invariant to the 
case of driven oscillators: 

E{t) = M [a(t)at(t) + at(t)a(t)] (3.1) 



d , , 



(See also [12], [IB], [3S] and [SS].) Here, A (t) = exp ^— J (c (s) — 2d (s)) ds^ and the correspond- 
ing time-dependent annihilation a (t) and creation (t) operators are explicitly given by 

1 p — 2ax ~ 6\ 
ait) = —iRx + e + t^ j, (3.2) 



aHt) = —\^px + e-t^ J (3.3) 

with p = i~^d/dx in terms of solutions of the Ermarov-type system fl2.5p - fl2.10p . These operators 
satisfy the canonical commutation relation: 

a{t)a^ (t) -a^ {t)a{t) = 1. (3.4) 
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The oscillator-type spectrum of the dynamical invariant E can be obtained in a standard way by 
using the Heisenberg-Weyl algebra of the rasing and lowering operators (a "second quantization" 
[28j, the Fock states): 

a (t) (a;, t) = (x, t) , (t) ^„ (x, t) = y/n + l (x, t) , (3.5) 

E (t) (x, t) = \ (t) (^n + ^„ (x, t) . (3.6) 
The corresponding orthogonal time-dependent eigenf unctions are given by 

'^n{x,t) = - , , ^ Hn{/3x + e), (^'^, ^'„) = (3.7) 

(provided that /3 (0) fi (0) = 1, when = A [21]) in terms of Hermite polynomials [39] and 

V'„(x,t) = e^(2"+^)^«vl>„(a:,t) (3.8) 
is the relation to the wave functions fl2.14p with 

^„(t) = -(2n + l)7(t) (3.9) 

being the Lewis phase [18], [26], [28] . 

The dynamic invariant operator derivative identity [6], [16]: 

can be veryfied in the following fashion. Introducing new linear momentum and coordinate operators 
in the form 

p=^(p-2ax-6), Q = X{l3x + e), (3.11) 

when [Q, P] = iX^ (a generalized canonical transformation), one can derive the simple differentiation 
rules 

^ = -2coaP'Q, ^ = 2aP'P. (3.12) 

(It is worth noting that if Cq = 0, the operator P becomes the linear invariant of Dodonov, Malkin, 
Manko and Trifonov [10], [11], [32], [55] for generalized driven harmonic oscillators.) 



Then ^ 

E = ^{P' + coQ') (Co = 0,1) (3.13) 

and it is useful to realise that E is just the original Hamiltonian H after the canonical transformation 
[26] . The required operator indentity fl3.10p can be formally derived with the aid of product rule 
(3.7) of Ref. ^6] (quantum calculus): 

and by fl3:T2|) : 

dE 

A— = Coa/32 (-gp -PQ + PQ + QP) = 0, (3.15) 
at 

which completes the proof. 
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Remark 1. The kernel 

is a particular solution of the Schrddinger equation li2.S\) for any solution of the Riccati-type system 
I[2.5\) - (2J1\) with Co = [4]. A direct calculation shows that this kernel is an eigenfunction 

{p - 2ax -6)K (x, y, t) = yK (x, y, t) (3.17) 

of the linear dynamical invariant 



4. Evaluation of Berry's Phase 

The holonomic effect in quantum mechanics known as Berry's phase [1], [2] had received con- 
siderable attention over the years (see, for example, |3], [12], [19], [20], [21], [22], [26], [35], [36] . 
[35] , [10], [H], [52], [5l] and references therein). The solution of the time- dependent Schrodinger 
equation (12. 2p has the form ( 12.16P with the oscillator-type wave functions (x, t) presented by 
(EH [21]: 

V^„(x,t) = e-^^"W vl/„(a;,t), (4.1) 

where ip^ (t) is the Lewis (or dynamical) phase and (x, t) is the eigenfunction of quadratic 
invariant (13.61) . (In the self-adjoint case, one chooses c = 2d when A = 1.) 

Then 

and Berry's phase On is given by 

^.idOn ^ A /■ ^* ^ ^ \ ^ A /^^^ ■ (4.3) 



rft \ Jr '' dt J ^"V \ -^t ", , 

Here, the eigenfunction is a 7- free part [26] of the wave function f l2.14p . namely 

and is, essentially, the real- valued stationary orthonormal wave function for the simple harmonic 
oscillator with respect to the new variable ^ = f3x + e (see (13.71) and (14.51) ). The integral (14.31) can 
be evaluated as in Refs. [22j and [261: 



, , ^ 5\E'„\ . / ^ f da n d6 dK\ ^ \ 1 , ^ / ^ d^n 

.da , .dS ,^ ^ , .dK ,^ ^ , 1 , / ^ (9$r 



i— ($„, x^^n) + i— x<l'„) + i— ($„, $„,) + - (c - 2rf) + ( $„, 



dt ^ ' dt dt 2 \ at 

where the last term is zero due to the normalization condition 



dx = 1. (4.5) 



00 



Moreover, 



($„,x^$„> = r'J^{e + e')^ld^ = r'(e' + n + ^^ 

/oo 
d^ = -er' 
-00 
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with the help of 



As a result, 



/ = 0, / di = n + 



OO -I 



dOn ^9/9 1\ da ^dS da 

dt V 2 dt dt dt 



(4.6) 



(4.7) 



and the phase On can be obtain by integrating (14.71) . Our observation reveals the connection 
of Berry's phase with the Ermakov-type system (I2.5I) - (I2.11I) . whose general solution is found in 
Ref. m. 



When c — 2d = f = g = 0, one may choose S = e = k, = and our expression (14. 7p simplifies to 



d6n 2 f l\ da 



(4.8) 



{n + - 

4a V 2 



,x2 cl' I ^ , f a! d' 



ji" jj, — iji) /i'/i + 2d\ /i 

a \ a d 



with the help of (12. lip . The function /i is a solution of the Ermakov equation (I2.12p -( l2.13p with 
Co = 1 and (3 = This result is consistent with Refs. [12] and [26], where the original expression 
of Ref. [36] has been corrected. 



5. An Alternative Derivation of Berry's Phase 



In view of and (gllD-dlJD, we get 

/ dOn d^n 
V dt dt 



Re {ij^, HtPn) = Re H^^) 



because the Hamiltonian in (I2.ip - (I2.2I does not involve time differentiation. Here, 



(5.1) 



H = ap^ + hx^ + - [px + xp) + - (c — 2d) — fx — gp 
and the position and linear momentum operators are given by 



X 



p 



/3 y/2a ,^ ^ 2ae 

^(a-at)+^(a + aT)+5-_ 



(5.2) 



(5.3) 



(5.4) 



in terms of the creation and annihilation operators (I3.2p - (l3.3p . After the substitution, the Hamil- 
tonian takes the form 

'a f 4a;^ 
2 



H 



h + 2ca i 

^ (c + Aaa) 

13^ 2 ^ [ 



(5.5) 



+ 



+ - 



a / 4a^ 
2 ['W 



^,1 b + 2ca i 



/3 

6 + 2ca 



[aa^ + a^a) + - (c — 2d) 
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+V2 



Aaa + c 



2ae 



^ (6 + ca) 



f + 2ga 

2/3 



+V2 



2 

Aaa + c 



+ i(^{g-2a6) + ^-{c + Aaa)] 



— I 



2/3 

73 



2ae\ 



—J-^(b + ca) 



2/3 



+a U - 



2ae 



+ 



/3 



/ + 



2a£: 



/3 



9 + 



ce 
1 



Here, 



- (a) 
2 ^ ' 



are the generators of a non-compact SU (1, 1) algebra: 

[Jo,J±] = ±J±, [J+,J- 



1 



Jo = - (aa^ + a^a) 



-2Jn 



(5.6) 



(5.7) 



and, therefore, a use can be made of the group properties of the corresponding discrete positive 
series I^-^ for further investigation of Berry's phase. (This is a 'standard procedure' for quadratic 
Hamiltonians — more details can be found in Refs. [3], [19], [29], [31], [33], [39], [H], [IE] and/or 
elsewhere.) Together, the linears and bilinears in a and a!^ realize the semi-direct sum of the SU (1,1) 
and the Heisenberg algebra fl3.4p (see Ref. [15] for more details). 



Thus 



ARe(^„,if^,) 



n + 



+a 5 



a ( /3^ + ^ 



+ 



h + 2ca 



(5.8) 



2a£ 



he 



2ae 



ce 



by (I331)-(I31D. 

Finally, from (13. 9p and (15. ip we arrive at a different formula for Berry's phase 

1~ 



~dt 



n + 



Aa^ h + 2ca 

- /3 + 



(5.9) 



2ae\ 

~ j ' /3 



+ 



/ + 



6£ 



5 



2ae\ 



9 + 



ce 



which is consistent with the previuous expression (14.71) for any solution of the Ermakov-type system 
(O-dnO]) (co = 1) . 
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